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1. Introduction 

Initially introduced in 1989 the Bose-Hubbard model has been the object of 
much recent work. It represents a simple lattice model of itinerant bosons which 
interact locally. This model turns out to describe fairly well recent experiments with 
bosonic atoms in optical lattices ^2jE]- Its low-temperature phase diagram has 
been uncovered in several studies, both analytical (see e.g. (QUEUED and numerical 
O ones. When parameters such as the chemical potential or the tunneling 
amplitude are varied the Bose-Hubbard model exhibits a phase transition from a 
Mott insulating phase to a Bose-Einstein condensate. Fig. |2j below, depicts its 
ground state phase diagram. 

In this paper, we investigate the phase diagram of this model in a mathematically 
rigorous way. We focus on the situation with a small tunneling amplitude, t, and a 
small chemical potential, [i. We construct the critical line between Mott and non- 
Mott behavior to lowest order in the ratio t/fi. More precisely, we prove the existence 
of domains with and without Mott insulator. These domains are separated by a 
comparatively thin stretch; the domain without Mott insulator is widely believed 
to be a Bose condensate. Our results establish in particular the occurrence of a 
"quantum phase transition" in the ground state. 

Over the years several analytical methods have been developed that are useful for 
the study of models such as the Bose-Hubbard model. They include a general theory 
of classical lattice systems with quantum perturbations (31 03 H3 ■ These methods 
can be used to establish the existence of Mott phases for small t; but they only apply 
to domains of parameters far from the transition lines. The Bose-Hubbard model on 
the complete graph can be studied rather explicitly and its phase diagram is similar 
to the one of the finite-dimensional model f|] . Results using reflection positivity are 
mentioned below and only apply to the hard-core model. A related model with an 
extra chessboard potential was studied in (see also [T7|). 
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The Bose-Hubbard model is defined as follows. Let A C "L d be a finite cube of 
volume |A|. We introduce the bosonic Fock space 

= H A ,N, (1.1) 

where Wa,tv is the Hilbert space of symmetric complex functions on A^. Creation 
and annihilation operators for a boson at site x € A are denoted by c\ and c x , 
respectively. The Hamiltonian of the Bose-Hubbard model is given by 

H = -t ctc y + ^U^2ctc x {clc x -l) (1.2) 

\x-y\=l 

The first term in the Hamiltonian represents the kinetic energy; the hopping param- 
eter t is chosen to be positive. The second term is an on-site interaction potential 
(assuming each particle interacts with all other particles at the same site) . The inter- 
action is proportional to the number of pairs of particles; the interaction parameter 
U is positive, and this corresponds to repulsive interactions. In our construction 
of the equilibrium state, we work in the grand-canonical ensemble. This amounts 
to adding a term — /j,N to the Hamiltonian, where N = c x c x is the number 
operator, and \x is the chemical potential. 

The limit U —* oo describes the hard-core Bose gas where each site can be occupied 
by at most one particle. This model is equivalent to the xy model with spin ^ in 
a magnetic field proportional to \i. Spontaneous magnetization in the spin model 
corresponds to Bose-Einstein condensation in the boson model. The presence of a 
Bose condensate has been rigorously established for fj, = (the line of hole-particle 
symmetry). See jj] for a proof valid at low temperatures in three dimensions, and ^3] 
for an analysis of the ground state in two dimensions. The proofs exploit reflection 
positivity and infrared bounds, a method that was originally introduced for the 
classical Heisenberg model in |llj . At present, there are no rigorous results about 
the presence of a condensate for \x ^ 0, or for finite U. 

The ground state phase diagram of the hard-core Bose gas is depicted in Fig. 
Q and reveals three regions: a phase with empty sites, a phase with Bose-Einstein 
condensation in dimension greater or equal to two, and a phase with full occupation. 
Particle-hole symmetry implies that the phase diagram is symmetric around the axis 
H = 0. 

The critical value of the hopping parameter in the ground state of the hard-core 
(he) Bose gas is 

too - ^. <u> 

This follows by observing that the cost of adding one particle in a state of vanishing 
density (where interactions are negligible) is —fi — 2dt. For fi < —2dt the empty 
configuration minimizes the energy, while for \i > —2dt a state with sufficiently low, 
but positive density has negative energy. The Mott phases of the hard-core Bose gas 
at zero temperature are stable because of the absence of 'quantum fluctuations' - 
the ground state is just the empty or the full configuration. The hard-core model is 
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Figure 1, Zero temperature phase diagram for the hard-core Bose 
gas. Bose-Einstein condensation is proved on the line fi = 0, for any 
t > 0. Our perturbation methods provide a quantitative description 
of the Mott insulator phases with density and 1. 



an excellent approximation to the general Bose-Hubbard model when t is small and 
/x is sufficiently small. 

A first insight into the ground state phase diagram of the general Bose-Hubbard 
model is obtained by restricting the Hamiltonian to low energy configurations. 
Namely, for — oo < \i < |£7, low energy states have or 1 particle per site. The 
restricted model is the hard-core Bose gas. Next, for < ll < |J7, states of lowest 
energy have 1 or 2 particles per site. The restricted model is again a hard-core 
Bose gas, but with effective hopping equal to 2t. We can define projections onto 
subspaces of low energy states for all fj,; corresponding restricted models yield the 
following approximation for the critical hopping parameter: 

.approx, , _/M if-00<//<itf, 

1 wkW) lf ( fc - v u < v < ( k + v u i 

(thin lines in Fig. |5J). The true critical line t c (fj>) agrees with ic Ppr ° x (/-t) up to 
corrections due to quantum fluctuations. We expect that 

t c (^) = t— M(l + 0(£)). (1.5) 

In order to state our first result, we recall that the pressure p(/3, fx) is defined by 

p03,M)= lim T^logTr^e-^"^. (1.6) 

K/% d |A| 

Here the limit is taken over a sequence of boxes of increasing size; standard argu- 
ments ensure its existence. Its derivative with respect to the chemical potential is 
the density; i.e., 

p(P^) = ~ P {P^)- (1-7) 
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Figure 2. Zero temperature phase diagram for the Bose-Hubbard 
model. Lobes are incompressible phases with integer densities. Thin 
lines represent the approximate critical line defined in p.4|) . 



The zero density phase is simpler to analyze because of the absence of quantum 
fluctuations. The following theorem holds uniformly in U, and therefore also applies 
to the hard-core model. 

Theorem 1.1 (Zero density phase). 

For n < —2dt, there exists (3q such that if (3 > (3$, we have that 

(a) the pressure is real analytic in [3, p; 

(b) p(/3,/i)<e-^. 

Here, a > depends on t, p, d, but it is uniform in (3 and U . 

This theorem is proven in Section |^1 

The transition lines between the Mott phases of density p ^ 1 and the Bose- 
Einstein condensate are much harder to study because of the presence of quantum 
fluctuations. We consider a simplified model with a generalized hard-core condition 
that prevents more than two bosons from occupying a given site. The Hamiltonian 
is still given by (jl.2j) . but it acts on the Hilbert space spanned by the configurations 
{0,1, 2} A . The phase dia gram of this model is depicted in Fig. This model is the 
simplest one exhibiting a phase with quantum fluctuations. Notice that, in the limit 
U — > oo, this model coincides with the usual hard-core model. The zero-density 
phase and the p = 2 phase are characterized by Theorem ll.il The transition line of 
the p = 1 phase is more complicated. The following theorem shows that it is equal 
to p/2d to first order in t/U, as in the hard-core model. 

Theorem 1.2 (Mott phase p = 1 in generalized hard-core model). 

Assume that < p < ^ and t < — const jj (with const ^ 2 11 d). Then there exist 

(3q and a > (depending on d, t, p,U) such that if (3 > /3q, we have that 
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Figure 3. Zero temperature phase diagram for the Bose-Hubbard 
model with the generalized hard-core condition. 

(a) the pressure is real analytic in (3, p; 

(b) <e"^. 

The critical line is expected to be close to t = for small i, /i, so that our 
condition agrees to first order in t. While we do not state and prove it explicitly, a 
similar claim holds around Adt = fi — U. Indeed, the p = 1 phase prevails for t < 
for small t. The "quantum Pirogov-Sinai theory" of references [3] applies here 
and allows to establish the existence of a Mott insulator for low t. Proving that the 
domain extends almost to the line t = requires additional arguments, however; 
Theorem ll.'il is proved in Section|31 The generalized hard-core condition considerably 
simplifies the proof. Indeed, it allows for a cute and convenient representation of 
the grand-canonical partition function in terms of a gas of non-overlapping oriented 
space-time loops, see Fig. |3]in Sectional The result is nevertheless expected to hold 
for the regular Bose-Hubbard model as well. 

While we cannot establish the presence of a Bose-Einstein condensate, we can 
prove the absence of Mott insulating phases away from the critical lines, by estab- 
lishing bounds on the density of the system. 

Theorem 1.3 (Absence of Mott phases). 

(a) For t > — ^ and for any A large enough, the density of the ground state is 
bounded below by a strictly positive constant, that depends on t, /i but not on 
A. This applies to the model with or without hard-core condition. 

2 

(b) Consider the model with generalized hard cores. 

Fort> 2d + Ct (h) d+2 and 
for any A large enough, the density of the ground state is less than a constant 

that is strictly less than 1; it depends on t, /i but not on A. 

This theorem is proved in Section It is shown that C ^ ^^(2 10 d-7r rf ) 3 +2 . 

Quantum fluctuations have some influence on the phase diagram, and a detailed 
discussion is necessary. "Quantum fluctuations" are fluctuations in the ground state 
around the constant configuration with k bosons at each site, for some k depending 
on /x and t. They are present in Mott phases for k ^ 1, while the ground state for 
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fx < —2dt is simply the empty configuration. Quantum fluctuations are not present 
in effective hard-core models where each site is allowed either k or k + 1 bosons. 
Their presence lowers the energy of both Mott and Bose condensate states. The 
key question is which phase benefits most from them. In other words, writing the 
critical hopping parameter as 

tc(M) = (M) (! + <#)), (1-8) 

the question is about the sign of a{jj), for small jj. 

The study in jlUj , based on expansion methods (no attempt at a rigorous control of 
convergence is made), suggests a rather surprising answer: the sign of a(jj) depends 
on the dimension! Namely, the quantum fluctuations favor Mott phases for d = 1, 
and they favor the Bose condensate for d ^ 2. We expect that this question can be 
rigorously settled by combining the partial diagonalization method of [S] with our 
expansions in Section |21 

2. Low-density expansions 

In this section, we present a Feynman-Kac expansion of the partition function 
adapted to the study of quantum states that are perturbations of the zero density 
phase. In this situation, quantum effects are reduced to a minimum, amounting 
basically to the combinatorics related to particle indistinguishability. Nevertheless, 
the resulting cluster expansion must deal with two difficult points: arbitrarily large 
numbers of bosons and closeness to the transition line. Both difficulties are resolved 
by estimating the entropy of space-time trajectories in a way inspired by Kennedy's 
study of the Heisenberg model ^B] — the present situation being actually simpler. 

The grand-canonical partition function of the Bose-Hubbard model is given by 

Z(P,A,n) =Tre- /3( ^ 7V) , (2.1) 

where the trace is taken over the bosonic Fock space. A standard Feynman-Kac 
expansion yields an expression for Z in terms of "space-time trajectories", i.e. 
continuous-time nearest-neighbor random walks. More precisely, 

N^O ' xi,...,x N eAneS N J J 

H exp{-C7 f 5g i{T)MT) dr}. (2.2) 

Here, 9 denotes a space-time trajectory, i.e. 9 is a map [0, 0\ — > A that is constant 
except for finitely many "jumps" at times < t\ < ■ • • < r m < (3, and 

\9(r j -)-9(r j +)\ = l. 

The "measure" Vxy on trajectories starting at x and ending at y introduced in Eq. 
(|2.2j) is a shortcut for the following operation. If / is a function on trajectories, then 

d^ y (9)f(9) = Y,t m £ / dn...dr m f{9). (2.3) 

m ^0 xi,...,x m -i J0<T 1 <---<T m <f3 

\xj—Xj-\\=\ 
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The second sum is over nearest-neighbor sites such that \x\ — x\ = |x m _i — y\ = 1. 
The trajectory 9 on the right side of 1)2. 3 j) is given by 

0(t) = x 3 - for j G [Tj,r j+ i), 

where (xq,t ) = (x,0) and (x m ,r m+1 ) = (y, /3). The underlying trace operation 
constraints the ensemble of trajectories to satisfy a periodicity condition in the "/in- 
direction". The initial and final particle configurations must be identical, modulo 
particle indistinguishablity. This explains the sum over permutations of N elements, 
7r G Sn, on the right side of (|2.2|) . 

We shall rewrite the expresssion (|2.2|) for the partition function in a form that 
fits into the framework of cluster expansions. The main result of cluster expansions 
is summarized in the appendix, and it is enough for our purpose. 

Trajectories are correlated because of (i) the interactions in the exponential factors 
of (|2.2|) which penalize intersections, and (ii) the permutations linking initial and 
final sites of different trajectories. The cluster expansion is designed to handle the 
former factors, but we need to deal first with the latter issue so to fall into the 
required framework. To this end, we concatenate each original trajectory with the 
one starting at its final site, so as to obtain a single closed trajectory that wraps 
several times around the (3 axis. Hence, instead of open trajectories [0, 0\ — > A, 
we consider ensembles of closed trajectories 0i : [0, £i/3] — > A, with ti being their 
winding number. Each such closed trajectory corresponds to a cycle of length ti of 
the permutation ir determined by the endpoints of the component open trajectories. 
For each cycle, the sum over ti sites in A and the integrals over the ti enchained open 
trajectories can be written as a sum over a single site Xj, followed by an integral 
over closed trajectories with 0j(O) = 9(£i(3) = Xj. Recalling that there are 

1 Nl 

permutations with k cycles of lengths £%,... ,1^, we obtain the following expansion 
of the partition function in terms of closed trajectories instead of particles: 

fc>0 X!,...,x k £A£i,...,e k ^l J 

k 

Y[w(6i) [J (l-CiKMi))- (2-4) 

Let 1(6) denote the winding number of the trajectory 6 : [0,£(9)/3] — > A. Its weight 
is defined by 

W {Q) = ^ye^W exp{ -UW (9) } . (2.5) 
Here, W(9) measures the self-intersection of 9, that is, 

w (°)= E / s e(ip+r),eup+r)dr. (2.6) 
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It will suffice to use the bound w(6) ^ e^rj e/3 ^ e "' '• Finally, interactions between 
trajectories 9 and 9' are given by 

Cu(0,9') = l-exp{-UW(9,9')}. (2.7) 
Here, W(9, 9') measures the overlap between trajectories 9 and 9', 

e(9)-l£(0')-l n 

W(e,9')=J2 E / $ew+T),e'W+T)dT. (2.8) 

i=0 j=0 Jo 

Expression (|2.4|) is suited for an application of Theorem IA.1I We show that the 
weights w(9) are small in the sense that they satisfy the "Kotecky-Preiss criterion" 

Proposition 2.1. For each closed trajectory 9 let j(9) denote the number of jumps 
of 9. Then, there exist constants a,b > such that 

/W r 

EVS dv2(e')e am+m 'tu(e,9>) < aj(9)+pbi(9). 

Proof. Since (u(9, 9') is increasing in U, it is enough to prove that, for any trajectory 
0, 

p f3(p+b)l' i- 

E E / ^ ^(0, 9') < aj(0) + W) - ( 2 -9) 

with 

^(0, 0') = X \e(iP + t) = 9'(k/3 + r) for some < i{9) - 1, 

O^fc^ £(0') - 1,0 < r < /3 

Here, X[-] denotes the characteristic function of the event in brackets. 

A trajectory 9' intersects 9 if a jump of 9' intersects a vertical line of 9, or if a 
jump of 9 intersects a vertical line of 9' (or both). Let denote the measure on 
trajectories [0, f3] — > Z d , starting at x = and with a jump at r = 0. Integration 
with respect to v§ can be defined similarly as in ()2.3j) : formally, we can also write 

j &4(9)f{9) = J dv$ (9)8(n)f(9), (2.10) 

where Vq is as in (|2.3jl (with x = y = 0), and where the Dirac function £(ti) 
forces the first jump to occur at n = 0. We get an upper bound by neglecting 
the restriction that trajectories need to remain in A. The left side of (|2.9|) is then 
bounded by 

/3(u+b)e' r , o 0(tJ.+b)e' r 

m E / <o* am + w) E !L V~ / d -o ^ am - (2.H) 

The first term accounts for trajectories 9' intersecting jumps of 9; the second term 
accounts for trajectories 9' involving a jump that intersects a vertical line of 9. We 
integrate over all trajectories [0, £'(3] — > 7L d that start at x = 0, without requiring 
them to stay in A. 
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Each trajectory 6' in the last integral of (|2,11|) can be decomposed into the jump 
from a neighbor z into 0, which contributes a factor te a (see definition ((231)), P 1 

us a 

trajectory from z' to 0. As has 2d neighbors we see that 

/^^>^/,? W( *. „ 12) 

Furthermore, the definition (j23J) implies that for every x, y, 
[ dvtf (6') ^ = £ (te a ) m E / dr i • • • dr - 

m^O xi,...,x m -i J0<Ti<-<T m <£'/3 

\xj-Xj-!\=l (2-13) 
\xi— x\=\y— a; m _i|=l 

< e te a 2d£'/3 

From (E3J), (gUJ) , (l2~T2l and (t2~T3l we conclude that 



< (j(9) + 2dte a f3£(6)) ^ exp{/3/[/x + 2die a + 6]}. 



(2.14) 



As /i + 2<it < 0, we can choose o and b such that fj, + 2dte a + b < 0. Then (|2.9f) holds 
for /3 large enough. □ 

Proof of Theorem \l.l\ Recall expression (|1.6|) for the pressure. Proposition 12.11 es- 
tablishes the convergence of cluster expansions, as stated in Theorem IA.11 With 
if denoting the usual combinatorial function of cluster expansions, see (|A.2|) . the 
partition function has the absolutely convergent expression 



a, = E E / d ^x(^) 

mjsl xi,...,a; m £A £i,...,£ m ^l 

m 

d^ m (^) ^ • • • m n ( 2 - 15 ) 



Taking the logarithm and dividing by the volume, standard arguments show that 
boundary terms vanish in the thermodynamic limit, and we obtain 

= E E E / d ^l) / d^J 2 (02) 

~ m 

■ J du^f x je m )< P (e 1 ,...,e m )l[w(e j ). (2.16) 



m>l x 2 ,...,x m gZ d fi,...,£ m >l ' 



Integrals can be viewed as functions of [3, fi, indexed by m, (xi), and (£«). They are 
real analytic in the domain (/3, /i) : (3 > (3o(fi). Their sum is absolutely convergent 
and Vitali's convergence theorem implies that p(fl, /x) is analytic. 
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Recall that the density is given by the derivative of the pressure with respect to 
the chemical potential; see (|1.7[) for the precise definition. The analyticity implied 
by the expansion allows for term-by-term differentiation. We can check that 



«/.) = Ef^'« ? fE"' E E 

„ m 

duif X2 (9 2 )...j dvi z ej6 m )< P (e 1 ,...,e rn )l[w(e j ). (2.17) 



Note that ^(0) = (3£(9)w(9), as follows from definition J23J) of the weight of 
trajectories. By (|A.5|) . we have the bound 

p((3,p) du^{9)w{9)e< 9 \ (2.18) 



There exists e > such that p + 2dte a + 6 + e < 0. Using l|2.13|) . we get 



L>1 



Then p((3, p) ^ e £ ^ for /3 large enough, and this completes the proof of Theorem 

rm ' □ 



3. Space-time loop representation 

The study of the transition line for the Mott phase with unit density requires 
the analysis of perturbations of the "vacuum" formed by one particle at each site. 
This involves the control of full-fledged quantum fluctuations. We turn, then, to 
a more general expansion setting previously employed to study spin and fermionic 
systems (HI • This setting shares some similarities with that of Section [21 but it 
also differs from it in significant ways. We use the same symbols u, w, £, £, 9, but 
we caution the reader that they are defined in slightly different ways. Besides the 
quantum-fluctuation issue, bosonic systems present the additional complication of 
the unboundedness of occupation numbers. In the present paper we wish to leave 
this second issue aside. We consider, thus, the model with generalized hard-core 
condition that ensures that configurations have at most two bosons at each site. 

Recall definition H2.1|) of the grand-canonical partition function. It is convenient 
to write 

H — pN = V + T, (3.1) 

where V denotes the diagonal terms (i.e., interactions and chemical potential terms) 
in the basis of occupation numbers in position space, and T denotes the hopping 
terms. We will consider T to be a perturbation of V. Our expansion is based on 
Duhamel's formula, 

e -P(V+T) = e -f3V + f' 3 dTe -rV { _ T)e -(P-rKV + T)^ ^ 
JO 
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which we can iterate to obtain 
c 



f3(v + T) = J2f dT l ...dr m ^ v {-T)^-^ v ...(-T)e-^) v . 

™^r, J0<ri<...<r m </3 

(3.3) 



Then 



Z(f3, A, n) = Tr e -^( y+T ) = ]T t m / dri . . . dr, 

m>0 JQ<Ti<-<T m <P 



• m 

'0<ri<...<r m </3 

£ TVe--^t iC , 1 e-(---^...4 m c J/m e-^--) y (3.4) 

(asi,j/i),...,(sDm,!/m) 

We denote by n = (n x ) x£ \, n x G N, a "classical configuration" that represents the 
state where n x bosons are located at site x, and \n) the corresponding normalized 
vector. Inserting projector decompositions 1 = £V. |nj)(rij| the trace can be written 
as 

Trp- r i y rt r P -(*a-n)V t r p -GS-rm)V _ 
J.I e L, Xl L. yi v . . . c Xm Ly m v — 

("o| e~ TlV '4i < Vi l n i)( n il e" {r2 ~ Tl)y . . . 4 m c ym |n m ) 

no,ni,...,n m 

As the operator V is diagonal in the base \n), this decomposition allows us to 
rewrite the expansion (|3,4j) in the form 

(/3,A,/x) = y d^(n) W (n) , (3.6) 

where 

(i) n is a space-time quantum configuration, namely an assignment of a config- 
uration n(r), for each < t < /3, such that 

— n is constant in r, except at finitely many times t\ < ... < r m , with m 
even. 

— At each tj, a "jump" occurs, i.e. there are nearest-neighbor sites (xi,yi) 
such that 

^x(Tj-) + l ifx = Xi, 
n r ( r ; + i - < n x (Ti-) - 1 if a: = y h (3.7) 
n x( T i— ) otherwise. 



Z 



(3.8) 



— n is periodic in the t direction: n(/3) = n(0). 

(ii) iu(n) are positive weights defined by 

w(n) = exp{ - f V{u{t)) dr} f[ [i^/n Xi (r i +)n ?/i (r i -) 
"'° i=l 

with the short-hand notation V(n) = (n| y |n). 

(iii) Integration with respect to the "measure" v on quantum configurations 
stands for a sum over configurations at time 0, a sum over m, integrals 
over jumping times, and sums over locations of jumps. 
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The expansion just obtained is rather general. It is convenient to interpret it in 
terms of random geometrical objects in a model-dependent fashion. For the case 
of interest here, we follow the "excitations", namely the sites where the occupa- 
tion number is different from the vacuum value 1. We therefore embed the "space 
time" A x [0, 0\ in the cylinder M. d x S 1 (with periodic boundary conditions in the 
time direction) and decompose the trajectories of the excitations in connected com- 
ponents. In this way, a quantum configuration n can be represented as a set of 
non-intersecting loops (with winding numbers n = 0, ±1, ±2, . . . ) in this cylinder. 
The representation is defined by the following rules: 

• The constant configuration n with n x (r) = 1, for all x € A and ^ r ^ /3, 
has no loops. 

• A jump of a boson from yi to Xi at time Tj (see (|3.7jl ) is represented by a 
horizontal arrow from (yi,Ti) to (xj,Tj). 

• The points (x, r) with n x (r) ^ 1 are represented by vertical segments. These 
segments point upwards if n x (r) = 2, and downwards if n x {r) = 0. 

Loops are illustrated in Fig. Similar representations have been used in various 
contexts, e.g. in a study of the Falicov-Kimball model JH]- Given a loop 7, we 



A 



Figure 4. Illustration for the gas of space-time loops. There are 
three loops with respective winding numbers 1,0, and -1. 



introduce the number of jumps ^(7) (always an even number, possibly zero); the 
length £0(7) °f all vertical segments pointing downwards; the length £2(7) of all 
vertical segments pointing upwards; £(7) = ^0(7) + ^2(7); and the winding number 
2(7). Notice that £2(7) — ^0(7) = z (l)- A loop 7 defines a unique quantum 
configuration n 7 . We define the weight of a loop as 

3(7) 

u,( 7 ) = tiM H JnKn-)nZ i (T i +) e - e °^e~ e2 ^ u -^ . (3.9) 

8=1 

Note that we have subtracted the classical energy of the background configuration 
with one boson at each site. The weight 10(7) thus only depends on excitation 
energies. 

These definitions allow us to rewrite the partition function (|3.6|) in terms of loops 
and their weights instead of space-time configurations. Unlike the trajectories of 
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Section the loops here have only a hard-core interaction due to the requirement of 
non-intersection. Furthermore, if T = { 7 i, . . . , 7m } is a set of disjoint loops, we have 
the important property that the weight of the corresponding quantum configuration 
n r factorizes, 



w (n r ) = e^l A l 



n«(7i)- (3-10) 

1=1 

We define a measure on loops, also denoted v, and we rewrite the partition function 
as 



1 



Z(l3,A,^) = e^J2— I Mil 



/m 
M-rm)Hw(<n) J] (i-C(7„7i))- 
i=l l<i<j<m 



ml , 

mJsO 1=1 l^i<j; 

(3.11) 

Here, the term corresponding to m = is set to e^' A ', and the function C(7>7') 
equals 1 if the loops 7 and 7' intersect (more precisely, if some of their vertical 
segments intersect), and equals if the loops have disjoint support. 

The expression (j3.11|) for the partition function is an adequate starting point for 
the method of cluster expansions. We prove that the weights are small so as to 
satisfy the "Kotecky-Preiss criterion", Eq. (|A.4[) . We can then appeal to Theorem 
lA.ll to conclude that the cluster expansion converges. 

Proposition 3.1. Under the hypotheses of Theorerri M.E. we have that, for any loop 

1, 

^(7X70^7, Y) e a(V) < a( 7 ) 



with 

a( 7 ) = 2 14 ^i(7) + 2 12 ^( 7 ). (3.12) 

Its proof relies on the bounds stated in the following lemma. Let us partition 
the set of loops into C = £( 0+ ) U where £( 0+ ) (resp. C^) is the set of loops 

with nonnegative (resp. negative) winding numbers. For each site z we introduce 
the measures v z on loops that make a jump at time r = involving z. Further, we 
let C z , C^ + \ and c[ ^ denote the sets of loops that contain (z,0). 

Lemma 3.2. Under the hypotheses of Theorem \l.°A for any site z, 

(a) J £(0+) d^( 7 )^(7)e a W < 2 n dg. 

(b) f 4 o +) di/( 7 ) w(rf) e a ^ < e -/3(C/-M-2 12 4) + 2 13^ 

(c) di^( 7 ) w( 7 ) e a W < Adte-^- 2dt - 2l2d2t ^\ 

(d) f,(-) Ml) w(rr) e a W < e -/3(M-2^-2 12 d 2 ^)^ 

Proof of Proposition ^ '.11 Suppose that the loops 7 and 7 ' intersect, i.e. £( 7 , 7 ') = 1. 
Then either a jump of 7 ' intersects a vertical line of 7 , or a jump of 7 intersects 
a vertical line of 7 ' (both may happen at the same time). The first situation is 
analyzed using the measures u z , and the second situation involves the sets df + ^ and 
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d \ More precisely, we have that 
di/( 7 , )w(7')C(7,7 , )e a(7 ' ) < 

(3.13) 

£(7) sup / d l / 2 (7 , )w( 7 ')e a ( 7 ') +j(7)sup / di^V) w( 7 ') e a{7 ' } . 



Using the estimates in Lemma 13.21 the right side is seen to be smaller than 0(7), 
provided [5 is large enough. □ 

Proof of Lemma \3. 6 A (a) and (b). Loops of £( 0+ ) have large energy cost, so crude 
entropy estimates are enough. Since ^2(7) ^ A) (7) for any loop 7 £ C^ + \ we have 
that ^£0(7) + (U- > Then 

^0(7) + {U- M )£ 2 ( 7 ) - 2 12 4£( 7 ) > \Ul{i). (3.14) 
Further, we can check that 

^Ht>/U> < 2 (3 15) 
From these observations and ()3.9|) . we obtain that 



i/7 



(7)e a{7) ^ < 



(4t)i(7) e -s w W if j > 2. 



(3.16) 



A loop with j(7) = 2n is characterized by a sequence of jump times ^ T\ < r 2 < 
• • • < T~2n- At each such time the trajectory can choose among at most 2d neighbors 
to jump to and 2 directions of time to proceed after the jump. The last jump is 
determined by the fact that 7 must be a loop, so there is no factor 2d (but both 
time directions are possible). The measure v z involves only loops with two jumps 
or more. From the last bound in (|3.16j) we obtain 

r / poo \ 2n— 1 

/ dv z (j)w(j)e a ^ sC Y 2-2 2n (2d) 2n ~ 1 (4t) 2n ( / dre"*^ 



(3.17) 



-1_(2£#)2- * 

Part (b) of the lemma follows from ()3.16j) and from considerations similar to 
(I3~T71) . Namely, 



40+) 



d^( 7 )u;( 7 )e a W ^ e-W-^d^) 



+ Y^2-2 2n (2d) 2n - 1 (At) 2n ( dre-i u A . (3.18) 



The first term in the right side represents loops without jumps. The right side is 
less than the upper bound in Lemma 13.21 (b^. □ 

Proof of Lemma \3. 6 A (c) and (d). Loops of have small energy cost when pa- 
rameters are close to the transition line. Estimates are needed that are more subtle 
than for loops of C^ + \ The situation is similar to that of Section [2 but a problem 
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needs to be solved: Loops, unlike trajectories, can backtrack in time. Our strategy 
is to first "renormalize" a loop 7 € by identifying a trajectory 6 = #(7) that 
moves only downwards, but with arbitrarily long jumps. Contributions of backtrack- 
ing can be controlled by similar estimates as above. The entropy of these trajectories 
can be expressed using an appropriate hopping matrix and we obtain sharp enough 
bounds. 

We start with (d). Given a loop 7 G Cq \ we start at (x,r) = (z,(3) and move 
downwards along 7. When reaching the end of a vertical segment (because of the 
presence of a nearest-neighbor jump), we ignore possible backtracking and directly 
jump to the next downwards vertical segment in the loop, at constant time. See the 
dotted lines in Fig. 0] We obtain a trajectory, since the motion is downwards only, 
punctuated by with long-range hoppings with which we must cope. 

Behind a hop from x to y there is a backtracking excursion between these sites. 
Its contribution to the total weight of the original loop (times e a ^ 7 ^) is given by the 
"hopping matrix" component 



a 



xy 



d^,(7M 7 )e a W, (3.19) 

where the integral is over loops that are open, have nonnegative winding number, 
start with a jump at (x, 0), and end at (y, 0). 

Each trajectory so constructed is characterized by a sequence of hopping times 
= T\ < ■ ■ ■ < T2m < P and a sequence of not-necessarily neighboring sites x = 
xq,xi,--- ,x m = x which are the successive hopping endpoints. Its weights are 
determined by factors exponentally decreasing with £q for each vertical segment and 
hopping matrix entries for each jump. In this way we obtain 

/ d,(7) -(7) e^> < e-^ 2 4) £ ^ £ f[ a 



— 1 

XQ=X 



m \-0 !W,-^m-lt=l (3.20) 



The overall exponential factor comes from the fact that £q ^ (3 because the winding 
number of the loops is not zero. The factor (5 m jm\ follows by integrating all choices 
of hopping times. 

To conclude, we must bound the sum of the matrix elements of a. The contri- 

2 14 d * 2 

bution of open loops that consist in just one jump is 2dte . Other open loops 
involve two jumps or more. Each jump has 2d possible directions. There are two 
possible time directions after each jump, except for the first and last ones. We need 
to integrate over time occurrence for each jump except the first one. We obtain 

2 / roc \ m—X 

Y^^ox < 2dte 2l4d ^ + ^2(2d) m 2 m ' 2 (At) m l / e^dr) 

x^O m>2 ' (3.21) 

< 2dte 2l * d £+2 9 d%. 
We used ()3.15|) . Inserting into (|3.2Uf) we obtain Lemma 13.21 (d). The bound of part 
(c) is similar, with an extra factor 2cfie u 1 ^ Adt for the additional first jump. □ 
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Proof of Theorem li.ifl This proof is similar to the one of Theorem 11.11 We use 
cluster expansions, in order to get a convergent expansion for the pressure, and prove 
analyticity by using Vitali's theorem. The density has an expansion reminiscent of 
(JXT7J), namely 

Jco dp 

di/(7 2 )... / di/(7 m )y>(7i,...,7 m ) JJw(7j). (3.22) 
The combinatorial function function <p is given by m> . From (l3~3]l 

^ = [^(7)-4(7)W7). (3.23) 

Again using Eq. (|A.5|) . we find the bound 

|p(/?,//)-l| < / dz.( 7 )|^( 7 )-4(7)h(7)e a(7) . (3.24) 

Only loops with nonzero winding number contribute. Going over the proof of Lemma 
13.21 (b) and (d) with 0(7) — > 0(7) + log £(7), we can check that the right side of the 
equation above is less than e~ £f3 whenever p — 2dt — 2 12 d 2 jj — e > and (3 is large 
enough. □ 

4. Density bounds 

Proof of Theorem M.cA (a). The Bose-Hubbard Hamiltonian preserves the total num- 
ber of particles, so that the density can be fixed. We denote by eo(p) the ground 
state energy per site in the subspace of density p. Neglecting repulsive interactions 
can only decrease the ground state energy; the minimum kinetic energy of a single 
boson is — 2dt. It follows that eo(p) ^ (— p, — 2dt)p for all £7^0. 

e(p) 




Figure 5. Upper and lower bounds for the ground state energy per 
site, eo(p). The density that minimizes e(p) necessarily satisfies p ^ 
a. 

We find an upper bound for eo(p) by using a variational argument. It is well- 
known that the symmetric ground state is also the absolute ground state, so that 
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we can consider a non-symmetric trial function. We decompose A into boxes of 
size I = L/0 1 / rf J . We consider the trial function 0f =1 (fj, where <pj is supported in 
the j-th box only and minimizes the kinetic energy. As is well-known, ipj is the 
ground state of the Dirichlet problem in the box, and the corresponding eigenvalue 
is — 2dt cos grj. Since £ + 1 ^ p~ x l d and cosx ^ 1 — this eigenvalue is less than 
-2dt + dt^p 1 ^) 2 . This implies that 

e (» ^ b{p) = (-/i - 2dt)p + 7r 2 dtp 1+ l (4.1) 

The minimum of b(p) is reached for p 2 l d = n z^^ 2 ) • m i n i mum value is 

'p + 2dt\ i+f 



By inspecting Fig. [5] we find that the ground state density is necessarily larger than 

» + 2dt ^ (4.3) 



d+2 \-K 2 t{d + 2 



□ 



Proof of Theorem \l.,% (b). The strategy is the same as for part (a), although quan- 
tum fluctuations bring extra complications. The variational argument leading to 
the upper bound for eo(p) can be modified by replacing particles with holes, so as 
to yield 

eo(p) < Hp) = -H+{n- 2di)(l - p) + Tr 2 dt{l - p) 1+1 d. (4.4) 

The lower bound is trickier. We fix the density p and work in the Hilbert space 
T~La,n with N = p\A\. We have that 

We can use the loop representation of Section |21 for the trace to obtain an expression 
similar to (|3.11|) ; the difference is that we require the sum of winding numbers of all 
loops to be equal to the negative of the number of holes M = |A| — N . 

The weights of loops with strictly positive winding numbers decays exponentially 
as e-W-ti, so they do not contribute in the limit (3 — > oo. We obtain an upper 
bound for Z(/3,A,p) (and therefore a lower bound for eo(p)) by neglecting the non- 
intersecting conditions between loops. Further, we replace the loops 7 with negative 
winding numbers by trajectories 9 as in the proof of Lemma 13.21 (c),(d). We then 
obtain the lower bound 

(4.6) 

Here, T denotes the multibody kinetic operator 

f = ^2,a{x-y)c\.Cy, (4.7) 

x,y 
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and a(x) is given in ()3.19j) . Then, by l|3.21|) . 

Jim |logTr ?iAM e-^ < [2dt + 2 10 d%)]M . (4.8) 

p — »-00 H 

The contribution of nonwinding loops is bounded using Lemma 13.21 (a), 

-i- f M-yM-r) < / ^o(lMl) < 2 u d£. (4.9) 
P|A| Jc(°) Jew 

We have shown that 

eo(p) >-n- (2dt -p + 2 10 d 2 £)(l - p) - 2 u d&. (4.10) 

From here on we proceed as before. The minimum of b(p) is — p— 7^2^572 { 2d ^~£" ) l+d / 2 . 
The ground state density then satisfies 

1- P > ^ d/Ad+2) l^L. (4.11) 

2dt-p + 2^d?^ 

One finds the condition of Theorem 11.31 by requiring that the numerator be strictly 
positive. □ 

Appendix A. Cluster expansions 

This appendix contains the main theorem of |2(J| for the convergence of cluster 
expansions. It allows for an uncountable set of "polymers", so that it applies here. 

Let (A, A, p) be a measure space with p a complex measure. We suppose that 
|/i|(A) < 00, where \p\ is the total variation (absolute value) of p. Let £ be a complex 
measurable symmetric function on A x A. Let Z be the partition function: 

Z = zZh. [ d ^--- [MAn) I] (l-C(^,A,)). (A.l) 

The term n = of the sum is understood to be 1. 

We denote by Q n the set of all (unoriented) graphs with n vertices, and C n C Q n 
the set of connected graphs of n vertices. We introduce the following combinatorial 
function on finite sequences {A\, . . . , A n ) of A: 

. s I 1 if re = 1 . . . 

tp(A 1 ,...,A n ) = \ (A.2) 

The product is over edges of G. A sequence (Ai, . . . ,A n ) is a cluster if the graph 
with n vertices and an edge between i and j whenever £(Ai,Aj) 7^ 0, is connected. 

Convergence of cluster expansion is garanteed provided the terms in (|A.1|) are 
small in a suitable sense. First, we assume that 

\l-C(A,A')\^l (A.3) 

for all A, A € A. Second, we need that the "Kotecky-Preiss criterion" holds true. 
Namely, we suppose that there exists a nonnegative function a on A such that for 
all A G A, 

J d\p\(A')\C(A,A')\e< A '^a(A), (A.4) 
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The cluster expansion allows to express the logarithm of the partition function as 
a sum (or an integral) over clusters. 

Theorem A.l (Cluster expansion). 

Assume that J d\/j,\(A)e a ^ < oo, and that (jA.3|) and ()A.4|) hold true. Then we have 
Z = exp|^ J d/i(Ai)... J dfi(A n )ip(A u ...,A n )y 

Combined sum and integrals converge absolutely. Furthermore, we have for all A\ € 
A 

1 + f d|/i|(A 2 )... I d\ f i\(A n )\i P (A 1 ,...,A n )\^e a ^\ (A.5) 

We refer to ^20 for the proof of this theorem, and for further statements about 
correlation functions. 
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